† Corresponding author. E-mail:
Project supported by the National Natural Science Foundation of China (Grant No. 81272495) and the Natural Science Foundation of Tianjin, China (Grant No. 16JC2DJC32200).
Recently, the phase compensation technique has allowed the ultrasound to propagate through the skull and focus into the brain. However, the temperature evolution during treatment is hard to control to achieve effective treatment and avoid over-high temperature. Proposed in this paper is a method to modulate the temperature distribution in the focal region. It superimposes two signals which focus on two preset different targets with a certain distance. Then the temperature distribution is modulated by changing triggering time delay and amplitudes of the two signals. The simulation model is established based on an 82-element transducer and computed tomography (CT) data of a volunteer’s head. A finite-difference time-domain (FDTD) method is used to calculate the temperature distributions. The results show that when the distances between the two targets respectively are 7.5–12.5 mm on the acoustic axis and 2.0–3.0 mm in the direction perpendicular to the acoustic axis, a focal region with a uniform temperature distribution (64–65 °C) can be created. Moreover, the volume of the focal region formed by one irradiation can be adjusted (26.8–266.7 mm3) along with the uniform temperature distribution. This method may ensure the safety and efficacy of HIFU brain tumor therapy.
In recent years, using ultrasound as a treatment modality has been increasingly widespread.[1–4] High-intensity focused ultrasound (HIFU) as a new technique to treat tumors has received a great deal of attention because it is noninvasive or minimally invasive. Also it can concentrate the ultrasound energy in a target area deeply seated in the human body and achieve the treatment purpose repeatedly. The HIFU has been used clinically to treat some solid tumors[5–7] such as uterine fibroid, breast carcinoma, prostatic cancer, etc. However, acoustic energy is difficult to deposit in the deeply seated brain tissue. The reason is attributed to the acoustic impendence of the skull. The absorption and attenuation of the skull are both strong. In early studies,[8,9] the researchers even removed part of the skull to enhance the ultrasound deposition on the target in the brain. With the development of the phased array transducer, ultrasound can propagate through the skull and focus into deep brain tissue using the adaptive focusing techniques.[10–14] However, its feasibility and safety are not ideal since temperature rise in the target area is difficult to control.
With the development of the phase compensation technique of sound wave and medical imaging technology, the HIFU is improved to penetrate the skull by researchers. In the 1990s, the research about the phase correcting technique showed that transcranial brain tumor treatment using HIFU was feasible.[11,14,15] From the end of the last century to recent years, the structure and parameter information of skull captured from medical imaging provided help for the simulation of transcranial HIFU.[12,16–18] The simulations and ex vivo verification experiments further demonstrated the feasibility of HIFU in transcranial brain tumor treatment. At the same time, studies showed that the hemispheric phased array transducer could maximize the penetration area of ultrasound on the skull surface and reduce the thermal deposition in the skull.[19] The proposition of hot spot eliminating algorithm and amplitude compensation helped to lower the temperature in the skull and raise the temperature at the target location.[10,13,20] The HIFU was implemented in some clinical trials to treat patients suffering from glioma, tremor, chronic neuralgia, etc. Some treatments were successful in partial tumor ablation,[21] reduced neuropathic pain,[22] and diminished tremor.[23,24] However, some treatments did not achieve the desired results.[22,24,25]
The focal temperature control in HIFU brain tumor treatment is critical. It affects the therapeutic efficacy directly. If the target area in the brain is overheated, it will be dangerous to the patient.[22,26] On the other hand, the treatment may be in vain if the focal temperature is insufficient.[23,27] In order to control the thermal deposition in the treatment volume, Ebbini and Cain used the pseudoinverse method to compute the array element amplitude and phase distributions to form multiple foci with desired levels at a set of control points in 1989.[28] Salomir et al. obtained a uniform temperature distribution within a large target volume using a double spiral trajectory of the transducer focal point in 2000.[29] Lu et al. proposed a multi-objective control method to generate an ideal spatial energy distribution for HIFU surgery.[30,31] Partanen et al. combined a multi foci sonication approach with a mild hyperthermia heating algorithm and achieved precise heating within the targeted region in 2013.[32] Zhou attempted to obtain uniformly distributed ultrasound energy in the target area by optimizing the irradiation path and time interval between irradiations in 2013.[33] The temperature distribution and volume of the focal region could be changed by electronic beam steering using a phased array transducer.[34,35] However, it is restricted more or less by the skull when HIFU is used to treat brain tumors.
In the present study, a new technique for modulating the transcranial temperature distribution of the focal region is proposed. In this technique, two targets are set to be either on the acoustic axis or perpendicular to the acoustic axis with a certain distance. As a result, two sets of excitation signals which focus on the two preset targets can be obtained using the time reversal method. The double signals are superimposed to stimulate each of the 82 elements of the transducer. The temperature distribution in the focal region is modulated by changing the triggering time delay and amplitudes of the two signals. Thus, a uniform thermal field is generated in the focal zone. The three-dimensional(3D) numerical simulation model is established based on an 82-element transducer and CT data of a 46 year-old male volunteer’s head. The simulations based on the FDTD method are implemented to calculate the thermal field formed by dual-signal superimposed triggering HIFU, also to test if the new method is able to build an ideal focal region with effective uniform temperature distribution.
Figure
In the simulations, the acoustic nonlinear propagation was described by the Westervelt equation written as[36,37]
The temperature distribution was calculated through the Pennes’ bio-heat conduction equation written as[38,39]
The equivalent thermal dose t43[40] was calculated from the following equation:
In this study, the parameters of skull and brain tissue such as ρ, c, α were obtained based on the high resolution CT of a 46 year-old male volunteer’s head provided by Tianjin Medical University Cancer Institute and Hospital. This study was approved by the ethics committee of Tianjin Medical University, Tianjin, China, and written informed consent was obtained from the volunteer. The CT scan parameters were 120 kV and 135 mA. The slice thickness and spacing were both 3 mm. To satisfy the spatial step of the simulation, a linear interpolation was performed between each slice of CT images.
The parameters were obtained from bone porosity (ϕ) converted from the Hounsfield unit (H) of the CT images and the calculation method was as follows:[16]
Here, F1 and F2 were set to be two different focal targets. The spacing between F1 and F2 was L as shown in Figs.
In this study, 55 °C is selected as the threshold of protein denaturation.[21,25] Unless otherwise indicated, the maximum temperature is no more than 65 °C since some side effects can be induced at over high temperature.[41] First of all, the simulations with the targets set to be on the acoustic axis (as Fig.
Taking the situation of L = 10.0 mm, Δt = 400 ns and M = 0.50 for example, the temperature distribution along the acoustic axis is marked as the blue curve shown in Fig.
The values of M for T1 = T2 with different values of Δt when the temperature reaches 65 °C are shown in Fig.
Figure
If ΔTm ⩽ 1 °C is selected as a necessity for the uniform temperature distribution, 0 < ΔzT ⩽ 6.5 mm can be obtained based on Eq. (
The modulations for generating a uniform temperature distribution are implemented for different L values, with the location of F1 fixed and the location of F2 varied. Figure
The modulation for generating a uniform temperature field in the direction perpendicular to the acoustic axis is implemented in this study. Take the case of F1 located at (0, 1, 75), F2 located at (0, −1, 75) and L = 2.0 mm (as shown in Fig.
Figure
The amplitude modulation is implemented based on the temperature distribution with Δt = 1000 ns as shown in Fig.
The targets are set as shown in Fig.
When the ultrasound wave propagates in skull, the mode conversion to shear waves is ignored in this study as most of the incident angles are less than 20° and the maximal incident angle is less than 25°.[10,43] The simulation model is established based on the CT data of a healthy volunteer’s head, which is justified by the fact that there is no significant difference in acoustic parameter between tumor and brain tissue.[44]
In this study, a method to modulate the temperature distribution and generate a focal region with uniform temperature distribution is developed. Two drive signals that focus on two preset different targets are superimposed. The trigging time delay and amplitudes of the two signals are adjusted to achieve the focal region modulation. Using a single focus of HIFU, the temperature distribution and volume of the focal region are able to be adjusted by changing the acoustic power or irradiation time.[45,46] On the other hand, it carries a potential risk of damaging the normal tissue on the path of ultrasound beam.[47–49] The temperature distribution and volume of the focal region can be changed in soft tissue by electronic beam steering or multi-foci sonication[28,32,34,35] using a phased array transducer. In this study, transcranial HIFU focus modulation is accomplished. The temperature distribution, length of long axis, short axis as well as the volume of the focal region are able to be modulated by this method. The uniform temperature distribution along the acoustic axis is realized by setting the minimal distance between the superficial target and the inner surface of the skull to be 10 mm. Also, the uniform temperature distribution on the y axis is achieved by setting the minimal distance between the superficial target and the inner surface of the skull to be 12.5 mm (Fig.
To improve the treatment efficacy of HIFU, Fan and Hynynen proposed that the temperature rise should be uniform in the focal region during the ultrasound irradiation.[50] Salomir et al. obtained a uniform temperature within a large target volume in the homogenous polyacrylamide gel and fresh meat samples using a double spiral trajectory of the transducer focal point.[29] Zhou attempted to obtain uniformly distributed ultrasound energy in the target area by optimizing the irradiation path and time interval between irradiations in his study using bovine liver.[33] However, the temperature modulation in transcranial HIFU has not been well studied. In this study, a uniform distributed temperature field was created in the brain by modulating delay time and amplitudes of two acoustic excitation signals. The focal region volume formed by one irradiation can be modulated effectively while avoiding overheating injury. The volume and temperature distribution of the focal region can be adjusted to meet the requirement of clinical treatment. Stacking the different sized focal region to cover the whole tumor is required in the treatment of a large tumor. The uniform temperature distribution along the direction of the acoustic axis can also be created when the targets were set both on one side of the acoustic axis and parallel with the acoustic axis (Fig.
The simulation results show that the range of L and Δ t for generating the uniform temperature distribution in the y direction are both narrow. The reason is that the length of the focal region in the vertical direction of the acoustic axis is short. The previous modulations of transcranial focusing were all completed based on the model as shown in Fig.
Connor and Hynynen indicated that no tissue exposed to less than an equivalent exposure of 30 min at 43 °C was damaged in the study on thermal deposition in the skull during transcranial HIFU treatment in 2004.[51] In Ding’s work about the modulation of transcranial focusing thermal deposition in HIFU brain surgery, the temperature rise at the skull was controlled to be less than 10 °C.[10] In this study, the temperature of the skull is controlled to be below 47 °C when the modulation is carried out to generate a uniform temperature distribution in the focal region. Theoretically, it will not cause skull and surrounding tissues to be damaged. However, the cooling water circulation system is essential during the treatment. Using the double excitation signal superimposition technique, the location of the focus is mobile, however, a sidelobe is not observed when the modulation is implemented to generate the uniform temperature distribution in the focal region based on the simulation model used in this study. Further study in our laboratory will focus on the consistency of this method. Also the probability and potential solution of the sidelobe will be investigated.
The HIFU transcranial therapy research has shown that the best way to enhance the therapeutic efficacy is to control the shape and volume of the focal region and meanwhile to modulate the temperature distribution in the focal region to generate an ideal uniform temperature distribution. In this study a modulation method is proposed by superimposing double acoustic excitation signals. A series of numerical simulations is performed to calculate the thermal field which is formed by HIFU through superimposing a dual excitation signal.
The simulated results are as follows.
i) The uniform temperature distribution can be created in the direction along the acoustic axis and in the direction perpendicular to the acoustic axis as well (Figs.
ii) The triggering delay time and the modulation coefficient of amplitude are two sensitive variables for modulating the thermal field (Figs.
iii) The distance between the locations of the peak temperatures is an important parameter to determine whether a uniform temperature distribution can be generated.
iv) When the distance between the two preset targets is 7.5–12.5 mm with the targets both on the acoustic axis, the focal region with uniform temperature distribution (64–65 °C, ΔTm ⩽ 1 °C) can be created by adjusting trigger time delay and the amplitudes of the two driving signals. The volume of treatable focal region formed by one single irradiation can be adjusted in a range of 26.8–95.3 mm3. When the distance between the two preset targets is 2.0–3.0 mm with the two targets on the opposite side of the acoustic axis, a focal region with a uniform temperature distribution can be generated by the same technique, and the volume of treatable focal region formed by one irradiation can be adjusted in a range of 88.0–226.7 mm3.
In summary, the numerical simulation results have verified that the new method proposed in this paper is capable of building a focal region with an effective treatable uniform temperature distribution. Also it can adjust the deposition volume of ultrasound energy during HIFU transcranial brain tumor treatment. We hope this method could improve the safety and efficacy of HIFU brain tumor therapy in the near future.
[1] | |
[2] | |
[3] | |
[4] | |
[5] | |
[6] | |
[7] | |
[8] | |
[9] | |
[10] | |
[11] | |
[12] | |
[13] | |
[14] | |
[15] | |
[16] | |
[17] | |
[18] | |
[19] | |
[20] | |
[21] | |
[22] | |
[23] | |
[24] | |
[25] | |
[26] | |
[27] | |
[28] | |
[29] | |
[30] | |
[31] | |
[32] | |
[33] | |
[34] | |
[35] | |
[36] | |
[37] | |
[38] | |
[39] | |
[40] | |
[41] | |
[42] | |
[43] | |
[44] | |
[45] | |
[46] | |
[47] | |
[48] | |
[49] | |
[50] | |
[51] |